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The class of the Generalized Coherent Potential Approximations (GCPA) to the Density Func- 
tional Theory (DFT) is introduced within the Multiple Scattering Theory formalism with the aim 
of dealing with, ordered or disordered, metallic alloys. All GCPA theories are based on a com- 
mon ansatz for the kinetic part of the Hohenberg-Kohn functional and each theory of the class is 
specified by an external model concerning the potential reconstruction. Most existing DFT imple- 
mentations of CPA based theories belong to the GCPA class. The analysis of the formal properties 
of the density functional defined by GCPA theories shows that it consists of marginally coupled 
local contributions. Furthermore it is shown that the GCPA functional does not depend on the 
details of the charge density and that it can be exactly rewritten as a function of the appropriate 
charge multipole moments to be associated with each lattice site. A general procedure based on the 
integration of the 'qV laws is described that allows for the explicit construction the same function. 
The coarse grained nature of the GCPA density functional implies a great deal of computational 
advantages and is connected with the 0{N) scalability of GCPA algorithms. Moreover, it is shown 
that a convenient truncated series expansion of the GCPA functional leads to the Charge Excess 
Functional (CEF) theory [E. Bruno, L. Zingales and Y. Wang, Phys. Rev. Lett. 91, 166401 (2003)] 
which here is offered in a generalized version that includes multipolar interactions. 

CEF and the GCPA numerical results are compared with status of art LAPW full-potential 
density functional calculations for 62, bcc- and fcc-based, ordered CuZn alloys, in all the range of 
concentrations. Two facts clearly emerge from these extensive tests. In first place, the discrepancies 
between GCPA and CEF results are always within the numerical accuracy of the calculations, both 
for the site charges and the total energies. In second place, the GCPA (or the CEF) is able to very 
carefully reproduce the LAPW site charges and a good agreement is obtained also about the total 
energies. 

PACS numbers: 31.15.Ew, 71.55.Ak, 71.23.-k, 71.15.Nc 



I. INTRODUCTION 

After forty years of studies and applications it is now 
clear that the density functional theory (DFT)^'^ con- 
stitutes a formidable tool for the understanding of the 
matter. Nowadays, DFT-based total energy calcula- 
tions^iii^!^ and Car-Parrinello molecular dynamics simu- 
lationsi are used in a growing number of scientific fields, 
ranging from physics to chemistry to biology. The rea- 
son of such an ubiquitous fortune is that these methods 
are ab initio, in the sense that the only underlying mod- 
els are the fundamental interactions laws and quantum 
mechanics. However, just because of their ab initio na- 
ture, DFT-based methods generally require large com- 
putational resources. In spite of the availability of faster 
and faster computers, this circumstance sets up the lim- 
itations to the applicability of the same methods. 

Most DFT implementations are based on the Kohn- 
Sham scheme^ and require the solution for the wave- 
functions of the appropriate Kohn-Sham Schroedinger 
equation. This usually implies the orthogonalization or 
the inversion of large matrices and, hence, a number of 
operations scaling, in principle, as N^, where N is the 
number of atoms in the system. While for semiconduc- 
tors or insulators, the wave- functions localization quite 
naturally leads to sparse problems, the case of metals ap- 
pears to be the most challenging. For metallic systems, in 



fact, the computational effort required by wave- functions 
based approaches remains 0{N^). Nevertheless, even 
for metals, approaches based on the direct minimization 
of the Hohenberg-Kohn functional wrt. the charge den- 
sity^i^i^ could achieve 0{N) scaling. In this case, the 
basic strategy consists in partitioning the system under 
consideration in a collection of weakly interacting frag- 
mentsi^ii^. Even with nowadays computers, the scaling 
properties of DFT algorithms wrt. the size of the system 
remain a crucial issue since they determine which classes 
of phenomena can be studied by ab initio methods. 

Among DFT implementations, the oldest methods us- 
ing such a 'divide and conquer' strategy for metallic sys- 
tems are perhaps the self-consistent versions of the Ko- 
rringaM, Kohn and Rostoker^- multiple scattering the- 
ory (MST). The MST methodS views the system under 
consideration as a collection of fragments (usually in a 
one to one correspondence to lattice sites) whose scatter- 
ing properties are determined by solving a Kohn-Sham 
Schroedinger equation. Once the fragment (or single- 
site) scattering matrices are determined, they are assem- 
bled together with the free electron propagator in or- 
der to obtain the scattering matrix, or, equivalently, the 
Green's function of the system. Both the determination 
of the fragment scattering matrices and the potential re- 
construction are 0{N), while, in principle, the solution 
for the global scattering matrix is an 0{N^) problem. 



2 



as it corresponds to the determination of the appropri- 
ate boundary conditions for the wave- functions in each 
fragment. However, a number of algorithms have been 
devisedi^ii^iii^ that are able to obtain 0{N) scaling by 
mapping the determination of the system's Green func- 
tion in a sparse problem. This is usually obtained by 
assuming zero the electronic propagator outside the so 
called Local Interaction Zone (LIZ) of each fragment. If 
the free electrons propagator is used, about ten neigh- 
bors shells should be included in the LIZ^^, while us- 
ing screened propagators^^ allows to have much smaller 
LIZ's: typically one or two neighbors shellsi^ are suffi- 
cient. Another remarkable feature of the MST method 
is that, being based on Green functions rather than on 
wave- functions, it can deal easily with disordered sys- 
tems and ensemble statistical averages. For this reason, 
since many years, the Coherent Potential Approximation 
(CPA) theory^° for disordered alloys has been used in 
conjuction with the MST^i and the DFTSS.. 

The present paper shall be concerned with the study 
of metallic alloys in which the nuclei are assumed to oc- 
cupy the positions of an ordered lattice, while substitu- 
tional disorder may be permitted. For these systems, in 
spite of the apparent complexity of the DFT algorithmic 
implementations, the analysis of large supercell calcula- 
tions has allowed for the identification of remarkably sim- 
ple trends. Namely, the charge excesses associated with 
each lattice site appear to be linear functions of the elec- 
trostatic potentials at the same site^^'^''. These simple 
relationships, to be referred in the following to as to the 
'qV laws, allow to describe the 'atoms' of each chemical 
species in an extended metallic system in terms of two 
parameters, say, the slope and the intercept of the above 
linear function o^^'^^ , and appear to be the appropriate 
generalization of Pauling's concept of electronegativity^! 
to solid state physics. We have already suggested that the 
'qV laws can lead to important simplifications for total 
energy calculations in metallic alloys^^. In the present 
paper, we shall introduce the class of the Generalized 
CPA's (GCPA) for dealing both with ordered and disor- 
dered metallic alloys. From the computational point of 
view, GCPA schemes present 0{N) scaling. Their prin- 
cipal virtue, however, is that, as we shall demonstrate, 
the GCPA functional exactly reduces to a function of 
the relevant charge multipole moments at the various 
lattice sites, thus constituting a coarse grained approx- 
imate version of the original DFT. At a further level of 
approximation, the GCPA density functional leads to a 
Ginzburg-Landau functional, the Charge Excesses Func- 
tional (CEF)^^, which is equivalent to the above linear 
'qV laws and computationally inexpensive. The predic- 
tions of the GCPA and the CEF about the 'qV laws 
and total energies shall be compared vs. full-potential 
Linearized Augmented Plane Waves (LAPW) calcula- 
tions^!^ for 62 ordered crystal structures^i^. Our con- 
clusions shall be that, at least for the systems considered, 
both GCPA and CEF are generally able to find out cor- 
rectly the system ground state and to fairly well repro- 



duce the energy differences between ordered structures 
in a fixed concentration ensemble. 

The following of this paper is organized as follows. In 
Sect. II we shall briefly review the MST version of the 
DFT. In order to have a functional form as much localized 
as possible, the relevant electrostratic contributions shall 
be rewritten using an exact multipole expansion. In Sect. 
Ill, we shall introduce the class of the GCPA theories and 
investigate the analytical properties of the corresponding 
approximate density functional. Moreover we shall ob- 
tain, as a further approximation to the GCPA, the CEF 
theory, already obtained in a much more phenomenolog- 
ical context^, that here is offered in a generalized form 
suitable for the inclusion of dipole or quadrupole inter- 
actions. In Sect. IV we shall compare the numerical 
results obtained from the GCPA and CEF approxima- 
tions with those from full-potential LAPW calculations. 
CEF and GCPA calculations appear numerically indis- 
tinguishable one from the other and both theories appear 
able to fairly well reproduce the LAPW total energies. In 
the final Sect. V, we shall draw our conclusions, make our 
comments and briefly discuss the possible developments 
of CEF and GCPA theories. 



II. REVIEW OF THE DENSITY FUNCTIONAL 
MULTIPLE SCATTERING THEORY 

A. The MST formalism 

In this subsection we shall briefly overview the grand 
canonical ensemble formulation of the MST-DFT— i^. Our 
aim shall be developing a common ground for dealing 
both with ordered and substitutionally disordered sys- 
tems. Although finite temperature, relativistic and mag- 
netic generalizations could straightforwardly be carried 
out^, in this paper we focus on the non-relativistic, non 
spin-polarized case at T = 0. Furthermore, when not 
otherwise stated, we shall consider the Local Density Ap- 
proximation (LDA)& to the DFT and assume to have ions 
of charge +eZi fixed at the lattice positions R^. 

In our discussion, the relevant density functional is the 
electronic grand potentia l^^'^^ . 
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where V is the volume of the system, ^ is the chemical 
potential and Etot is the sum of the total electronic 
energy and the nuclei electrostatic interaction. N{£] /i) 
is the integrated density of states which is related to to 
the electronic density of states (DOS), n(e,/i), through 
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the following relationship: 

N{e, fJ-) — I de n{e, fx) 



(2) 



The notation highlights the implicit /i dependence of the 
DOS that arises from the effective Kohn-Sham potential. 
In a frozen ions treatment, of course, the nuclear inter- 
actions term is just a constant that is included here for 
future convenience. 

The basic idea underlying the MST is partitioning the 
system in 'small' scattering volumes, Vi, X^i = 
which in most implementations are 'centred' at the nu- 
clei positions. Although at this stage the partitioning is 
quite arbitrary, as we shall see in the following, there is a 
natural choice for it. Using the Lloyd's formul a'^^i'^^ , the 
integrated DOS, iV(e; /i), can be expressed as the excess 
with respect to the corresponding free electrons quantity, 
iVO(e): 



N{e]^i) = N°{e) - -ImlndetM{e) = 

TT 

7V°(e)-f-/mVTr(lnT(£)),, (3) 

TT ^ — ' — 

i 

where the trace is taken only over the angular momentum 
components. In Eq. ([3]) the multiple scattering matrix, 
M , or the scattering-path matrix'^^, r — M_~^ are defined 
in terms of the single-site scattering matrices^, t^ (e) , and 
the free electron propagator, G^j{s), is given by: 
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It is convenient to recall here that the single-site scat- 
tering matrices convey the informations about the phase 
shifts at the surfaces delimiting each scattering volume. 
The continuity of the wave- functions at the same surfaces 
is ensured by the construction of the scattering-path ma- 
trix, T, this is accomplished by the numerical inversion of 
the multiple scattering matrix M. Since the size of M is 
proportional to the number of scatterers in the problem, 
its inversion is the source of 0{N^) scaling in the MST 
version of the DFT. 

Within MST the link between the electronic density 
and the scattering matrices is provided by the Green 
function^ 

Gij.LL'{r,r' ,e) = Z.,^L{r,e)Tij,LL'{£)Zj,L'{r',e) 
{r~r')Z,^L{r,e)Jj.L'{r',e) 

+eir' - r)ALir,e)Z,,L'ir',e)\SLL'S,j (5) 

where r e Vi, r' e Vj. ^^^^(r,^) and Ji^L{j^,£) are, re- 
spectively, the regular and irregular at r = solutions^^ 
of the KS Schroedinger equation for the energy e. For 
real energies both Zi^L{r,e) and Ji^i(r,e) are real func- 
tions. The (site resolved) charge densities, the DOS and 
the neat charges at the i-th site can be obtained by in- 
tegrating the Green function over the energy and/or the 



appropriate volumes and by taking the trace over the 
angular momentum indexes as follows: 

Pi{i^;fJ-) = / de^Im{G^^^LL{r,r' ^ r;e)} (6) 

ni{e;n) ^ [ dry^ Im{Gu.LL{r,r' ^ r;e)} (7) 



As it is shown in Refs.l36l.l22l.l32l. the Hohenberg-Kohn 
density functional, Eq. ([T]), can be more conveniently 
rewritten within the MST formalism as the sum of a ki- 
netic and a potential energy functionals, as follows: 



n{T = 0,V,ii) = T - fxN + U 



(8) 



where the above two contributions are given by the fol- 
lowing expressions: 
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The effective potential in Eq. v'^-^^{r), is specified 
by the Kohn-Sham equation. 



dt 



v^^ir, [p]) 



'I ^Ir-R, 



(11) 



It consists of the Coulombian potential due to the elec- 
tronic and ionic charges and of the exchange-correlation 
potential, v^'^{r,[p]) = dE^^[p]/Sp{r), where E^^[p] 
is the third term on the RHS of Eq. dTU]). In the lo- 
cal density approximation (LDA) v^"-^ is assumed to de- 
pend locally on the electronic density, i.e., ^'^'-^(r, [p]) = 
v^'^{p{r)) and e^'^(r, [p]) = e^'^{p{r)). 

We wish to highlight a useful consequence of the above 
partitioning of the system volume. The density func- 
tional defined by Eq. ([8]), which is, of course, variational 
wrt the global charge density, p(r), turns out to be vari- 
ational also wrt the charge densities in each scattering 
volume Uj, in formulae. 



sn 



= 



Furthermore, it is possible to show'^'^ that 



SU 



(12) 



(13) 
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and that 



S{T - fiN) 



_ e// 



(14) 



It is interesting to observe that, the expression for the 
site resolved DOS, Eq. ([7]), aUows to recast the integrated 
DOS and the electronic grand potential as sums of site 
resolved contributions. These contributions, however, in- 
volve the site-diagonal part of the system Green function 
or scattering matrix, Gjj, or r^j, and then are non triv- 
ially coupled together through the boundary conditions. 
If this coupling was neglected, as it is done, for instance, 
in the case of the Harris-Foulkes density functionali^, 
0{N) scaling could be obtained. Fortunately, similar nu- 
merical performances can be achieved with less dramatic 
approximations. A sensible alternative is to impose ran- 
dom boundary conditions at the fragments surfaces^^. In 
this paper we shall follow a different approach and use 
averaged boundary conditions. As we shall see in the fol- 
lowing Section, this allows to have a tractable form for 
the coupling in the kinetic part of the density functional 
and permit to obtain 0{N) algorithms. Although it was 
proposed with a different aim, one of the oldest method 
applying such mean boundary conditions is the CPA, a 
generalized version of which shall be offered in the next 
Section. Before, however, we need to discuss a differ- 
ent source of coupling that is present in the potential 
energy part of the functional, namely, the electrostatic 
interactions between fragments. In the past this subject 
has received little consideration and it has been ruled 
out by invoking the screening properties of metals. How- 
ever, nowadays there is a general consensus that careful 
estimates of these interactions are necessary in order to 
obtain accurate total energies for metallic alloys. 



B. 



Multipole expansions for the effective potentials 
and the potential energy 



We have shown in the previous Sect. Ill Al that the DF- 
MST theory is variational wrt the local charge densities, 
Pi(r), of each fragment or scattering volume. In this sub- 
section we shall see how the multipole expansion used by 
most numerical implementations of the theory has the 
conceptual advantage of giving expressions for the effec- 
tive Kohn-Sham potentials and the potential energy in 
which different scattering volumes are coupled together 
only through simple functions of the multipole moments. 

The relevant formulae can be obtained by splitting the 



volume integrals in Eqs. ([TO)) and (jlip in sums of inte- 
grals extending over the scattering volumes, Vi and by 
expanding the denominators in spherical harmonics. Al- 
though they require some labour, the derivations are very 
straightforward and need not to be reported here. The 
resulting expressions for the potential energy and the ef- 
fective potentials are listed below: 



MAD 
L 



(15) 



and 



«i (r) 



^2yMAD 



\r — r' 
(r) 



(16) 



In Eqs. p5|) and p6)) we have introduced the local 
multipole moments, 



qi,L = / drpL{r)pt{r) - ZiSl^^q^ 



0) 



and the Madelung potentials. 



MAD 



where 



^rMAD 
Vi J. 



j^i L' 



(17) 



(18) 



(19) 



The coefficients Mijxu = ^Ihu^R-ji) are given by 



Mii,(R) = 47r ^ 

L" (l"=l+l') 



E 



c 



^„ (2r + i)!!yi.(R) 



(2r-fi) 



(20) 

C^i^, are the Gaunt numbers^ and the functions pl{y) 
in Eq. are defined as. 



47r 



(2^ -hi)!! 



(21) 



The only values that are relevant for spherical approxi- 
mations are poo(r) = 1 and Moo,oo(^) = 1/-R- 

In Eq. (Uni), the contribution from the i-ih. lattice site 
to the potential energy is denoted as Mi([pi(r)]) and given 
by 



w»([a(i")]) = y 



dr I dr 



,Pi(r)p,(r') 



dr 



e^Zjp,(r) 



drp,(r)e^^(p,(r)) 



(22) 



Within the LDA, Ui depends on the electronic density 



at the i-th fragment only, while for non-local approxima- 
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FIG. 1: Triangular inequalities that must be satisfied in order 
to have a convergent multipole expansion: r < |r — (R2 — 
Ri)|, r < IR2 — Ri|. The partition of the system volume in 
Voronoi polyhedra, marked by the lines, guarantees that the 
inequalities hold. 



tions to the DFT there could be some dependence on the 
density at sites j ^ i. 

Much published work has been done within spherical 
approximations (SA), namely the muffin-tin (MT) or the 
atomic sphere approximation (ASA). In that case only 
the first terms, ^ = 0, of the multipole expansions are 
included. Thus, Eqs. p5|) and ([16]) must be replaced by 
the following expressions. 



U 



E 



MAD 



and 



(r) = 



dv 



+ t;^^(p,(r)) + e^y, 



(23) 



MAD 



(24) 

Thus, for SA's, the only relevant multipole moments 
are the local charge excesses. 



% = 9i,oo = / dTpi(T)-Zi 

J Vi 



(25) 



and the Madelung potentials are constant within each 
scattering volume to the values 



yrMAD _ T/AfAD 



^ R 



(26) 



Remarkably, in Eqs. (fT51) and or in their SA coun- 
terparts, Eqs. ([23]) and (j24|) . the charge densities at dif- 
ferent sites, Pi{r), are coupled only with the Madelung 
potentials at the same sites, V/^^^-'^{r). Of course, the 
last quantities contain information about the charge den- 
sities at all crystal sites. 

We wish to highlight that the multipole expansion does 
not converge for arbitrary partitions of the system. Actu- 
ally, convergence requires that, for any pair of scattering 
centers, and Rj, and for any point r belonging to the 
scattering volume Vi, the triangular inequality illustrated 



in Fig. (1) must be satisfied. It is easy to realize that par- 
titioning of the system in Voronoi polyhedra accordingly 
with the Wigner-Seitz construction guarantees the above 
condition to be fuUfilled everywhere except but for the 
zero measure set of points constituted by the surfaces 
of the polyhedra, thus ensuring the convergence of the 
theory. The Wigner-Seitz construction, therefore, con- 
stitutes a natural choice for the partitioning. 



III. GENERALIZED COHERENT POTENTIAL 
APPROXIMATIONS (GCPA) AND CHARGE 
EXCESSES FUNCTIONAL THEORY (CEF) 

A. Generalized Coherent Potential 
Approximations (GCPA) for the scattering matrices 

In this Section we shall discuss a whole class of approxi- 
mations for systems with atoms lying on a regular lattice, 
where, however, substitutional disorder is allowed for. 
Metallic alloys, both ordered intermetallic compounds 
and random alloys, constitute the most relevant exam- 
ple of such systems. Other examples are crystals with 
empty, or 'vacancy', sites. Although, in general, these 
systems do not have translational invariance, neverthe- 
less the underlying 'geometrical lattice' does. Forty years 
ago, this consideration led Soven to formulate the Coher- 
ent Potential Approximation or CPA^". Since then, the 
CPA had an appreciable fortune. Its crucial virtue was 
that, by introducing a 'mean field' fashion effective crys- 
tal, it allows to use many techniques designed for ordered 
systems that were already well developed at the time at 
which the theory was proposed. 

For many years, the DFT implementations of the 
Qp^22,38 j^ave been based on the assumption (in the fol- 
lowing referred to as the single-site approximation or SS) 
that sites occupied by atoms of the same chemical species 
are characterized by the same effective Kohn-Sham po- 
tentials. Although the DFT-SS-CPA has been proved 
able to carefully determine the electronic structure and 
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, nev- 



the spectral properties of many alloy systems 
ertheless it leads to an incorrect description of the elec- 
trostatics and of the total energies in metallic alloys^. 
Due to its mean field nature, in fact, the SS approxima- 
tion neglects the fluctuations of the charge transfers and 
the energetic electrostatic contributions associated with 
them. This failure has stimulated many authors that en- 
visaged CPA generalizations aimed to include the effects 
of different chemical environments^2ii2iiMiiSiii&. 

In this paper we define a class of approximations for 
DFT-based electronic theories in which most of the above 
CPA generalizations can be included. We shall refer to 
the approximations belonging to such a class to as Gen- 
eralized CPA (GCPA). A theory belonging to the GCPA 
class shall be identified by: (a) a theory specific 'external 
model\ i.e. a rule for determining the effective Kohn- 
Sham 'site' potentials and the statistical weights Wi to 
be assigned to each 'site', and (b) an approximate form 



for the kinetic part of the density functional, specified by 
Eqs. (|27ll30p below. The last feature is common to all the 
theories belonging to the GCPA class. 

Before discussing the ansatz for the kinetic functional 
we wish to illustrate what a GCPA 'external model' can 
be on the basis of a few examples. The first example of 
a GCPA theory is, of course, the DFT implementation 
of the SS-CPA in Refs. Q2k and [s^- Its external model 
is the SS assumption (identical effective potentials for 
atoms of the same atomic species and weights propor- 
tional to the respective atomic concentrations). Another 
example is the Polymorphous CPA (PCPA) of Ujfalussy 
et alj^i^i^. The external model is constructed using 
an auxiliary supercell containing TV atoms, usually hun- 
dreds or thousands, each to be weighted with the same 
weight. The effective site potentials are reconstructed 
on the same supercell via Eq. (jlip. thus atoms of the 
same chemical species are allowed to have different po- 
tentials depending on their environments. This specific 
choice for the external model appears the reason why 
the PGPA theory substantially improves the alloy elec- 
trostatics while maintaining all the advantages of the 
standard SS-CPA about the spectral properties'^. Other 
existing CPA-based approaches like, e.g., the Non-Local 
GPA^O'Si, ^j^g siM-CPAiii^ can also be considered as 
particular cases of GCPA's. 

We shall now introduce the kinetic ansatz that is com- 
mon to all GCPA theories. For this purpose we prefer not 
to start from the definition of the functional. Rather we 
shall follow a path closer to physical intuition and to the 
spirit of Seven's original CPA formulation^^. At similar- 
ity of SS-CPA calculations, the GCPA defines an effective 
periodic crystal whose sites are occupied by effective 'co- 
herent' scatterers characterized by the single-site scatter- 
ing matrix i'^(e), the corresponding Green function shall 
be Then, if we considers the Green function of 

a single substitutional impurity with a single-site scat- 
tering matrix embedded in the above effective crystal, 
Gjj(<j,i'^), the GCPA consists in requiring that 



(27) 



In other words, the weighted average of the impurity 
Green functions must be equal to the 'coherent' Green 
function G'^(t'^). In Eq. l[?f|) the energy dependences have 
been dropped for sake of simplicity and iV stands for the 
number of different scatterers in the model. 

Eq. ([77)1 is illustrated in Fig. (2). In terms of the 'co- 
herent' scattering-path matrix of the effective lattice, r'^, 
and of the CPA 'projectors', 1?^, it can be rearranged as 
follows: 



TV—* - 



(28) 
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FIG. 2: A pictorial illustration of the SS-CPA for a binary, 
Ac^Bcg, alloy (top) and of the GCPA (bottom). The rectan- 
gular frames have the meaning "the Green function of what 
is inside". 



A GCPA theory is then an approximation for the t 
matrix, whose diagonal elements are given by 



(30) 



while the diagonal matrix elements of the Green function 
are given by Eq. ([5]) with i = j. 

Within the approximation defined by Eqs. (|27ll30p 
above, the MST reviewed in the previous Section al- 
lows to calculate the charge densities and the integrated 
DOS, N{e,fj,) and, through Eq. ([5]), the kinetic part of 
the Hohenberg-Kohn functional. Here we need not to 
trace all the intermediate steps that can be reproduced 
following the scheme of Refi^. The GCPA approximate 
version of the Lloyd formula, Eq. is given by: 

= — — 1 — Im y WiirmT[e) 



N 



—Im 



(31) 



It has the very remarkable property that the integrated 
DOS N{e; n) and, hence, the kinetic functional are vari- 
ational^ wrt. both If and r"^. 

In Sec. II A, we have mentioned that in the exact MST 
the contributions to the integrated DOS associated with 
each lattice site and proportional to {Inr^n, are coupled 
together because each element of the r matrix depends on 
the scattering properties of all the lattice sites. Within 
the GCPA, the only source of coupling is t^. Each lo- 
cal contribution depends on and on the local poten- 
tial. However, within the GCPA, the Lloyd formula does 
not depend on nor on the local potentials. As a con- 
sequence, the integrated DOS results in a sum of local 
contributions, coupled together only through if , 



(32) 



We shall call this very controlled and tractable kind of 

(29) 

coupling marginal coupling. 
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In view of further developments, it is convenient to iso- 
late in Eq. ([5^ two distinct terms. The first arises from 
the first two addends in Eq. ([3T|) . it is identical for all 
sites and related to the effective background defined by 
the GCPA medium. The second depends on the local 
CPA projectors and, through them, on the local poten- 
tials. In formulae: 

N,{e■,^l) = iV^^'=(e;/i) + -ImTr\ogD^{e) (33) 

TT 

Implementing the GCPA within the DFT gives for the 
kinetic functional of Eq. ^ the following marginally cou- 
pled form: 



where 



and 



(35) 



T^i[p^],^i) = -ImTr de\ogD,{e) 

drpi{r;fi)vf^{r;^i) (36) 



As mentioned in Sect. II, in MST-based DFT calcula- 
tions, the only source for 0{N^) scaling is the inversion 
of the multiple scattering matrix, Eq. ^ , required to ob- 
tain the scattering-path matrix r. This step is bypassed 
in a GCPA theory by approximating the relevant matrix 
elements r^j via Eq. ([28]) in terms of the local scattering 
properties and the coherent scattering matrix jf, the last 
of which is, in turn, obtained by an averaging process. 
For this reason, GCPA theories are 0{N), allowing for 
very substantial savings of computing time. Of course, 
the price for these savings is payed by the approximation 
implied by Eq. (|30p . A diagrammatic analysis of these 
errors can be found in Ref. la- 
It is necessary to make a couple of remarks about the 
physical meaning of the GCPA in the present context and 
to highlight the differences with respect to the traditional 
way in which CPA-based theories have been introduced 
in the past. In first place, the GCPA has been introduced 
here as an approximation for the Hohenberg-Kohn den- 
sity functional. As an approximation, it may well be used 
to describe an ordered alloy. Its range of applicability is 
by no means confined to the realm of random alloys. In 
second place, the introduction of the weights Wi to be 
assigned to each scatterer makes GCPA theories suitable 
for dealing with sophisticated pictures of the order (or 
the disorder) in metallic alloys. This, of course, requires 
what we have called an 'external model'. In a foregoing 
paper we shall discuss a (to some extent) self-consistent 
way to define an external model that is able to provide 
a picture of ordering phenomena in metallic alloys as a 
function of the temperature. 



B. The DFT-MST-GCPA functional: the 'marginal 
coupling' property 

In the present subsection we shall analyze certain for- 
mal properties of the GCPA approximations introduced 
in the previous subsection. All the above discussion 
can be summarized in the following approximate density 
functional: 



GCPA 



Wi 



GCPA 



MAD 
i,L 



(37) 



In Eq. ([STjl the qi^L are defined by Eq. pT]) and the local 
part of the GCPA functional by 

c.f^^^([p,],Ai) = mp,],fi) + u,{[p,],fi) (38) 

where the terms Ti{[pi], fi) and Ui{[pi\, p) are given by 
Eqs. JSSl) and ^ above. We note that the local GCPA 
functional, wf^*-^^"^, depends also on the atomic number 
of the atom at Hi, Zi, and on the volume and the shape of 
the i — th Voronoi polyhedron through the local potential 
energy term, Ui. In the following we shall make the sim- 
plifying assumpion of having identical Voronoi polyhedra 
for all the sites considered. 

In Eq. ([37|) the couphng potentials, Vj^'[^^ , are pro- 
vided by the specific external model. In the following of 
this Section we shall assume: 



(39) 



Appropriate choices of the coefficients Al and of the 
weights, Wi, give then the SS-CPA or the PCPA. Fur- 
thermore, Eq. ([55]) can also be used for spherical or for 
full-potential charge reconstructions. 

As mentioned in Sect. II A, Eq. (fT2|) . the density func- 
tional is variational not only wrt the global charge den- 
sity, p(r), and the chemical potential p, but also wrt the 
charge densities in each scattering volume, Pi(r). More- 
over, as discussed in Sec. Ill A and in Ref. [32, the GCPA 
density functional is variational wrt the effective medium 
scattering matrix, r"^. Furthermore, in a GCPA theory, 
the background kinetic term, T'""=^{p) in Eq. §7} de- 
pends on the electronic density only through rf and p. 
Thus, the functional derivation of Eq. ([57]) wrt the local 
densities, Pi{r), gives the following set of coupled equa- 
tions: 



GCPA 



MAD 



(r)=0 



(40) 



where we have used Eqs. (HH) and (EH). 

Within a GCPA theory, solving the set of the Euler- 
Lagrange equations PO)) . one for each scattering cen- 
ter, together with the equations that determine the 
chemical potential and the coherent scattering matrix 
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r'^, is completely equivalent to the minimization of the 
density functional. As it is apparent, these Euler- 
Lagrange equations are coupled each other only through 
the Madelung potentials, and /i. Moreover, the func- 
tionals Luf^'^'^'^{pi{r)) are identical for sites occupied by 
the same chemical species. 

In order to understand the consequences of the above 
result let us consider, for instance, an alloy sample con- 
stituted by a large supercell. One may wish to calculate, 
in the given sample, the properties of different 'atoms', 
in first place the charge densities, Pi(r). Inside the sam- 
ple, ]f, and the cell geometry are fixed, thus the set 
of the Pi{r) is completely determined by the values of 
the Madelung potentials y/^'^^(r) and by the atomic 
number of the ion at the position R^. More gener- 
ally, inside the given sample, any site diagonal property 
Hi shall be completely determined by Zi and by the set 
of Madelung potentials. We can establish this result as 
follows: 



n, = n(z„\/,^'^^^(r)) (41) 

Examples of such site diagonal properties are the local 
contributions to the grand potential, the multipole mo- 
ments, and the local DOS. The functional forms, one for 
each alloying species, given by Eqs. ([1T|) sometimes can 
be easily numerically fitted and then constitute a useful 
tool for the evaluation of the site quantities in the given 
sample. They are the source of the simple laws, as e.g. 
the 'qV laws, empirically found from extended metallic 
systems calculations. This notwithstanding, GCPA the- 
ories are able to predict complex trends for certain site 
diagonal properties as, e.g., the site resolved DOS's^i^. 

Since Eqs. (|¥T|) allows to evaluate, among other prop- 
erties, also the charge density of each fragment and, 
hence, the full charge density, p(r), then, in virtue of 
the Hohenberg and Kohn theorem, it follows that any 
ground state observable in the sample given is a func- 
tional of the set of the Madelung potentials at all the 
crystal sites, Vj'^'^^^{r), and of the set of the atomic 
numbers only. Since the last is, again, specified by the 
sample, it follows the theorem: any ground state observ- 
able in the sample given is a functional of the Madelung 
potentials only, or, equivalently, a function of the set 
of coefficients, {V^^"^^}, that completely determine the 
Madelung potentials. 

Since, in virtue of Eq. (|39p . the coefhcients in the set 
I y MAD I g^j,g linear functions of the set of the multipole 
moments, {q}, then the above theorem implies the corol- 
lary that any ground state property of the sample is a 
function of the same moments. Within the GCPA and 
for the specific sample given, it is then possible to refor- 
mulate the DFT in terms of the charge multipole mo- 
ments. By neglecting a constant term with the physical 
meaning of the grand potential contribution due to the 



mean GCPA 'atom', Eq. ((37|) can be written as 

i 

e^ 

-i^^mqi,oo (42) 

i 

In deriving Eq. (|42|) . we have used Eq. (f39|) and the 
fact that, since wf'*-^^"^ is completely determined by the 
local density, Pi(r), it cannot depend on the multipole 
moments at other sites. Moreover, we have isolated the 
contribution proportional to the chemical potential, /x. 
Having introduced an explicit dependence on /x, the last 
term in Eq. (|^^ can be thought as a way of enforcing 
the global electroneutrality. This is unnecessary if we 
consider a specified sample in which, of course, /i has a 
precise, fixed vale. However, since the term proportional 
to fj, has precisely the form it must have, its introduc- 
tion is equivalent to extend the validity of Eq. (j42|) to 
all samples specified by the same mean atomic concentra- 
tions and by the same value for jf. 

To summarize, we have established the following re- 
sults. Within the GCPA class of approximations the 
Hohenberg and Kohn density functional can be recast 
in the form of Eq. (HH). It consists of (a) local terms, 
-GCPA -j-Qj, -j-j-^g scattering site, consisting in func- 
tions of the charge multipole moments that are identical 
for sites with the same chemical occupation; (b) a bi- 
linear form coupling the charge multipole moments at 
different sites, with coupHng coefficients Mij^LL' defined 
by the crystal geometry; (c) a term proportional to the 
chemical potential that ensures the global electroneutral- 
ity. The functional defined by Eq. (|^^ is identical for 
all the alloy samples characterized by the same mean 
atomic concentrations and the same value for the coher- 
ent scattering-path matrix jf. Evidently, it constitutes 
a coarse grained version of the DTP because the math- 
ematical definition of the multipole moments, Eq. (|17|). 
does not completely determine the charge density. The 
last is determined by the multipole moments only within 
the GCPA theory. This reduction of the relevant infor- 
mation has been obtained at the price (a) of the GCPA 
approximation and (b) of having restricted the considera- 
tion to a specific sample. Nevertheless, no restriction has 
been made about the size of the sample that, therefore 
can be chosen in such a way to guarantee an appropri- 
ate description for a fixed concentration ensemble, as we 
shall discuss at the end of the present section. 

Having recast the GCPA functional as a sum of func- 
tions of the charge multipole moments has obvious math- 
ematical advantages. However, we have not yet com- 
pletely determined the functional form of the local ener- 
getic contributions uj^'~'^'^{{qi}, Zi). In order to do this, 
we need to make the hypothesys that, in the sample con- 
sidered, the distribution of the Madelung potentials coef- 
ficients, {V^^^^}, is continuous in the range of the values 



that the same potentials assume in the sample. This is 
consistent with the observations in Refs. [23l. \24. [26l. Let 
us consider two scattering sites, say i and j, occupied 
by the same chemical species, a, at which the Madelung 
coefficients take very close numerical values, Vj^l^^ = 
Vl'"^^ and y/^^^(r) = yA^^^(r) + Ay*^^^(r). The lo- 
cal energetic contributions, the charge densities and the 
local multipole moments shall be: uuf''-'^^ = lu^'^^^, 
Pi{r) = p(r) and qi l = 9l for the i-th. site, and 

^GCPA ^ ^GCPA ^ ACj^CPA^ ^^.(j.) ^ ^(j.) ^ ^^(j.) 

'lj,L = Ql + Agi for the j-th site. To the first order in 
A/c»(r) we have: 



(t^ rt^) ft^ (ic) 




Alu 



-GCPA 



dry*^^^(r)Ap(r 



Ap(r) 



where Eq. (|40p has been used. The substitutions of the 
expansion for the Madelung potential, Eq. (fT5|) . and of the 
expressions for the charge multipole moments, Ea. (jl7p . 
then give: 



a<d: 



-GCPA 



MAD 



Aqi 



(43) 



Once integrated over qL Eq. (|43l) gives 



tc) (tc) ftc) (tc^ ) ( tc 



FIG. 3: The Local Interaction Zones (LIZ) used in PCPA 
(left frame) and LSGF (right frame) calculations are marked 
by the dark areas. For each scattering site, the kinetic func- 
tional is evaluated using the appropriate single-site scattering 
matrices, ti, inside the LIZ, and the effective CPA scattering 
matrix, t'^, outside. 



the purpose of the present paper, we notice that most 
non-local density functionals offered in the literature^ 
are actually local wrt the density gradients, thus most 
non-local schemes will remain marginally coupled in the 
above sense. 

(ii) Splitting the kinetic functional into local contribu- 
tions marginally coupled trough the coherent scattering 
matrix jf is a simplification due to the GCPA. In fact, it 
has been obtained by assuming averaged boundary con- 
ditions at the surfaces of the Voronoi polyhedra through 
Eq. ([77]) . An estimate of the so induced errors can be 
obtained by the comparison of PCPA vs. Locally Self- 



r,GGPA 



-GCPA 



({9°}) 



L 



/ V^^^^ {{q/})dq' Consistent Green Function (LSGF) calculationsii exe- 
J{qO} cuted on the same supercell. As it is sketched in Fig. (3), 



Eq. 

'qV 



(44) 

(|44|) can be easily numerically evaluated from the 
data, y/i^^ = Vl^^°{{q}) obtained as an out- 
put from GCPA calculations. Unless a constant with the 
meaning of the local energy contribution at {q} = {q^}, 
it determines the local energies for each chemical species 
a. Eqs. ((43|) and (j44|) have been obtained under very 
broad conditions: the differentiability of the kinetic func- 
tional^^'^'' and the monotonicity of the 'qV laws. The 
first is the usual requirement for the convergence of the 
Kohn-Sham scheme of the DFT, while the second con- 
dition is certainly verified by all GCPA calculations re- 
ported in the literature, including those executed at ex- 
tremely high values for the Madelung potential (see Fig. 
7 in the next subsection and the related discussion). 

In the remainder of the present Section we shall make 
a few general comments about the validity of the frame- 
work defined by the GCPA theory in comparison with 
the exact density functional. 

(i) The fact that the effective potential and the po- 
tential energy functional can be decomposed in site con- 
tributions coupled together only through the Madelung 
potentials is an exact consequence of the LDA and it 
has nothing to do with the GCPA. This kind of coupling 
is 'marginal' in the sense that, although not necessarily 
small, it has the simple and tractable functional form 
which arises from the bilinear terms involving the mul- 
tipole moments in Eq. (|42p. Although this is beyond 



cuted on the same supercell. As it is sketched in Fig. (3), 
both calculations evaluate the kinetic contribution from 
the i-th to the functional by solving the problem of a sin- 
gle impurity, in the case of the PCPA, or of an impurity 
cluster, the Local Interaction Zone (LIZ), for the LSGF. 
In both cases the scattering matrices outside the LIZ are 
set to the coherent scattering matrix. If. PCPA calcula- 
tions can then be viewed as LSGF calculations with only 
one atom in the LIZ. This argument also suggest that, 
wrt GCPA calculations, exact DFT results include, for 
each site, corrections depending on its chemical environ- 
ment. 

(iii) We have already seen that the coarse grained ver- 
sion of the GCPA functional, Eq. ([42]), holds for all the 
alloy configurations characterized by a specified value for 
t'^ in a fixed concentration ensemble. This could appear 
as a serious limitation since it looks unlikely that, e.g., 
T*^ could have the same functional energy dependence for 
two different systems. In general, in a GCPA theory, 
t'^ is a ground state property determined not only by 
the mean concentrations but also by the distributions of 
the Madelung potentials for each alloying species. As 
opposite to the SS model where these distributions are 
trivial, more sophisticated external models, as e.g. the 
PCPA, give complicated charge and Madelung potential 
distributions. How could then the GCPA functional be 
useful in such cases? As argued by Faulkner et ali^, 
the PCPA theory applied to ideal random alloys gives 
well-defined values for all physical properties. This is be- 
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FIG. 4; The PCPA supercells shown in frame (a) and (b) 
contains, respectively, A*' — > oo atoms in a random alloy con- 
figuration and n atoms in an ordered configuration, both at 
the same mean atomic concentrations. N is large enough to to 
guarantee an appropriate description of a random alloy within 
a GCPA theory. Similarly, n has been chosen to permit the 
description of an ordered alloy up to some length scale The 
supercell in the frame (c) is identical to that in (a) except 
but for the dashed region which contains n atoms in the same 
ordered configuration as in (b). The cell (c) is therefore able 
to describe an ordering fluctuation up to the scale set by I. 



cause, in perfect random alloys, the distribution of the 
chemical environments is easily obtained by statistical 
considerations and it is given by the appropriate multi- 
nomial distributions. Therefore, the PCPA random alloy 
constitutes a privileged reference system whose physical 
properties, including jf, can be approximated up to an 
arbitrary accuracy by letting the number of atoms in the 
PCPA supercell, N, going to infinity (see Fig. 4a). We 
believe that the same t'^ obtained for a random alloy at a 
given concentration can be used for building a physically 
clean, though approximate, theory also for ordered alloys 
at the same concentration. In the next section we shall 
provide numerical evidences for that, here we present a 
more formal argument. Imagine that an ordered array 
containing n atoms (Fig. 4b) is able to account for the 
properties of same ordered alloy configuration, up to a 
length scale, I, that can be made large at will in the 
n — > oo limit. In Fig. 4c we draw a supercell, a part of 
which is constituted by the supercell of Fig. 4b, while the 
remaining N — n sites are occupied as in the random alloy 
supercell of Fig. 4a. We can think that the supercell in 
Fig. 4c represents a fluctuation of an ordered phase in 
a random alloy matrix and that it describes the physical 
properties of such fluctuation up to the same length scale 
I as in Fig. 4b. We are implicitly using the common idea 
of 'locality' in physics, or, in a more specific context, of 
'nearsightedness' of the DFT— . However, as Eq. (|28|) im- 
plies, the difference between the coherent scattering-path 
matrices corresponding to Figs. 4a and 4c, — r^, is 
proportional to the ratio n/N and, then, it can be made 
small at will in the iV — > cxd limit, for any value of n. 
We conclude that the coherent scattering-path matrices 
of a random alloy, r^, is able to account for the physi- 
cal properties of ordered configuration considered. The 
limitation n/N << 1, that comes from the above argu- 
ment, does not impose any upper bound on the maximum 
length scale at which chemical fiuctuations can be stud- 
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FIG. 5: 'qV relationships for a bcc random Cuo.50Zno.50 
alloy. The excesses of electrons, qi = g,;,oo are plotted vs. the 
Madelung potentials, Vi^^'^ . The results have been obtained 
from Imax = PCPA calculations for a supercell containing 
432 atoms at lattice costant a — 5.50 a.u.. Circles represent 
Cu atoms and triangles Zn atoms. Note that positive values 
for Qi correspond to negative net charges and vice versa. 



ied and it is of no practical importance provided that N 
is large enough to ensure a good approximation r^. As 
reported in the literature*^'^^, it seems that N about 100 
is already enough. 

iv) Although we have suggested that the coherent scat- 
tering matrix from random alloys GCPA calculations can 
be used for ordered alloys too, we are aware of the limi- 
tations of such a physical picture. For instance, a GCPA 
theory always implies finite quasiparticle lifetimes^, and, 
hence, a smearing of the peaks of the Bloch spectral func- 
tion (BSF). 



C. A generalized version of the Charge Excess 
Functional Theory 

As a matter of fact, the analysis of DFT supercell cal- 
culations for metallic alloys suggests the existence of sim- 
ple relationships between the charge excesses at the lat- 
tice sites, Qi.oQ, and the Madelung potentials at the same 
sites, K- 



MAD 



Namely, simple linear laws, one for each 



alloying species, have been found to hold, say 



MAD 
i,00 



(45) 



where and ki have the same numerical values for atoms 
of the same chemical species in the given supercell. Ex- 
amples of the linear 'qV laws obtained from PCPA cal- 
culations for a binary and a quaternary alloy are reported 
in Figs. 5 and 6. 

It is interesting to observe that similar linear 'qV 
laws can be derived starting from the GCPA functional. 
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Eq. p2| . for random alloys by a second order series ex- 
pansion about the zero Madelung field multipole mo- 
ments, {q'^}, that can be obtained by solving the fol- 
lowing set of equations: 



duj. 



GCPA 



dqi. 



(46) 



FIG. 6: 'qV relationships for an fee random 
Alo.25Cuo.25Nio,25Zno.25 alloy. The exeesses of elec- 
trons, Qi — gi.oo are plotted vs. the Madelung potentials, 
y^MAD rpj^g results have been obtained from £max ~ 
PCPA ealculations for a bee supercell eontaining 108 atoms 
at lattice costant a = 6.88 a.u.. Circles, squares, triangles 
and crosses stand for Cu, Ni, Zn and Al atoms. 



This procedure leads to a Ginzburg-Landau configura- 
tional 'Hamiltonian' in which the relevant fields arc con- 
stituted by the values of the multipole moments of each 
lattice site. In formulae: 



i,L,L' i,j.L.L' i 



where we have omitted the term n^^'' {{q°}, = 0) 
that represent the GCPA energy at zero Madelung field 
and chemical potential and that is constant in a fixed 
concentration ensemble. The coefficients ai^L^i are given 
by the second derivatives of the GCPA functional 



a-i.LL' 



a~GCPA 



{{<l^}) 



(48) 



The functional of Eq. ((47|) constitutes a generalization 
of the Charge Excess Function (CEF) proposed in Ref.|23 
for discussing the charge transfers in metallic alloys and 
shall then be referred in the following to as the CEF. 
The novel feature here is that Eq. (|T7| includes not only 
the charge excesses, (/i^ooj but also the charge multipole 
moments with £ > 0. ^ 

The minimization of the CEF functional ^'-^^^ wrt. 
its variables, the set of the multipole moments, {g}, and 
the chemical potential, /i, gives 

^ [ai,LL'{qiX' - <7°L') + MLL',ij qj,L'] = Af5j,oo (49) 
L' 

and 



^ qim = 



(50) 



Using the definition of the Madelung potentials, 
Eq. ([31]), and setting 



^0 
LL'qi.L' 



00 



(51) 



it is easy to show that Eq. ([^5]) for L = (0, 0) coin- 
cides with the linear laws given by Eq. (|45|) . Versions 
of Eqs. (P^ . (|5D|) and (|?T|) with the angular momentum 
summations truncated at £ = can be found in Ref. .261 
We wish to highlight that the CEF derivation from 
the GCPA functional is based on the assumption, com- 
mon to all Ginzburg-Landau theorie a^^i^^ , that the ho- 
mogeneously disordered phase, in the present case the 
random alloy phase, can the starting point for a pertur- 
bative treatment of ordering or segregation phenomena. 
As discussed in the previous subsection, in the GCPA 
context, this amounts to conjecture that the coherent 
scattering-path matrix jf of a random alloy can be used 
for obtaining a physical picture of concentration fiuctua- 
tions, or, in other words, that it is able to represent such 
fluctuations. 

We wish to close this Section with a few comments. 
The principal result of this paragraph, the CEF func- 
tional of Eq. (|T7|) , has been obtained by a series expan- 
sion of the GCPA functional about the values that the 
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FIG. 7: 'qV relationships for light substitutional impuri- 
ties (vacancies, H and Li atoms) dissolved in bcc Al, from 
CPA+LF calculations. The quantity Q is the number of va- 
lence electrons at the impurity site, therefore Q = Q corre- 
sponds to electrons for vacancies and H atoms and 2 elec- 
trons for Li. The linear behaviors observed for small fields 
{Vmad < 1 a.u.) are superseded by power law trends (see 
the log-log plot in the inset) for very high fields. 



multipole moments would have in absence of coupling. 
The series has been terminated at the lowest order at 
which differences with respect to SS approximations are 
expected for. This, not surprisingly, is enough to obtain a 
physical picture of the charge transfers in metallic alloys. 
For a given alloy configuration, the linear Euler-Lagrange 
equations obtained by minimizing the CEF can be easily 
solved for the charge multipole moments. The procedure 
require the inversion of the matrix F_ of elements 



F, 



ij.LL' 



0,1, LL' Sij + Mij^LL 



(52) 



As we have shown elsewere^'^*' , for a given alloy config- 
uration, the value of the CEF functional at its minimum 
has the physical meaning of the total energy of the same 
configuration. The ambiguity due to the presence of the 
above mentioned concentration dependent constant can 
be resolved by comparing CEF and GCPA calculations 
for a single configuration in a fixed concentration ensem- 
ble. 

In the previous subsection we have described a general 
procedure based on the numerical integration of the 'qV 
laws for evaluating the functional form of iijf"~'^'^{{qi}). 
Of course, if the random alloy t'^ was able to represent 
concentration fluctuations and the 'qV laws were lin- 
ear, the GCPA and the CEF functionals would be co- 
incident. We do not think that the 'qV laws can be 
truly linear. The argument is as follows. The local ex- 
cesses of electrons, Qi^o, accordingly with the physical 
intuition and with the results plotted in Figs. 5 and 6, 
are non-increasing functions of the Madelung potential 
V^m'^- If the 'qV laws were really linear, qi_oo would 



decrease indefinitively and eventually reach unphysical 
values, Qi^oo < Zi corresponding to negative charge den- 
sities. Actually we expect that the linear laws cannot be 
any longer valid when all valence electrons are expelled 
from the site. This circumstance would correspond to 
some critical value for the charge excess, say g.f qq. Before 
this critical value is reached, the 'qV laws should exhibit 
a crossover to an asymptotic behavior, say gi.oo — > 9f oo 
as V/^f^^ oo. We have tested this conjecture by exe- 
cuting CPA+LF calculations^^ for a single impurities, va- 
cancies, H or Li atoms, embedded in Al. The results are 
shown in Fig. 7, where we plot Q = qifiO—qZoo for the im- 
purity site as a function of the relevant Madelung fields. 
In all the cases considered, a linear regime is clearly vis- 
ible at low fields. A very high fields, a crossover to a 
power law dependence is observed, with the number of 
electrons tending to the critical value from above. The 
crossover field is comparable with the host band width. 
We recall that in CPA-I-LF calculations t'^ is that of the 
host while the Madelung potential is just an adjustable 
papameter. While this is a sensible way for studying the 
response of the impurity to the perturbing field, this do 
not imply that all the range of the perturbations con- 
sidered is physically meaningful. We do not think that 
such high fields, corresponding to the tunneling regime 
in the impurity site, could occur in real systems as this 
would require a too large defect of electrons at the impu- 
rity nearest neighbors. Hence, Fig. 7, while supporting 
the view that the linearity of the 'qV laws and the CEF 
are just approximations, does not support the possibility 
that, at least for metallic systems, appreciable deviations 
from linearity or failures of the CEF are likely to occur. 

Another point we wish to address is concerned with the 
value of the chemical potential fj, in Eq. (|T7)l . In a recent 
paper, Drchal et al.^^ argued that fi should be always 
zero since the Fourier transform of the Madelung coeffi- 
cients with L ~ L' = (0, 0) diverges as fc — > cx) implying 
that the sum of the charge excesses ^ Qi.oa must vanish, 
automatically satisfying the electroneutrality constraint. 
The observation of Drchal et al. is correct for infinite 
systems, while for finite supercells, even with periodic 
boundary condition, the same Fourier transform always 
remains finite, k, in fact can take only the values of the 
reciprocal space vectors that consitute the tiling of the 
supercell considered^. The set of the allowed values for 
k includes only for infinitely large supercells. In most 
practical calculations, then, fi is necessary, although usu- 
ally it takes small non-zero values. 



IV. NUMERICAL RESULTS 

In this Section we present a series of numerical tests 
designed to study the limits of validity of the GCPA 
and CEF theoretical frameworks. The central issues here 
shall be investigating the realm of validity of the linear 
'qV laws, Eq. (HS)) or and of the energetics implied 
by the CEF functional, Eq. (ITTl) . Furthermore, we shall 
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try to answer two questions: (i) to what an extent the 
CEF is able to approximate GCPA calculations and (ii) 
how do the predictions from the CEF and the GCPA 
compare vs. 'exact' DFT calculations for ordered sys- 
tems. The GCPA theory chosen for these tests is the 
PCPA*'^, that, being based on a supercell approach, al- 
lows for easy comparison vs. 'exact' DFT calculations. 

Several kinds of calculations shall be presented in this 
Section. The 'exact' DTF results used for comparison 
shall be LDA full-potential LAPW calculations produced 
using the WIEN2K ab initio packag o^^i^^ . They are re- 
ferred below to as LAPW. In all the cases about lO'' 
k-points in the full Brillouin zone have been used, the 
spherical harmonics expansion of the potentials in the 
mufhn-tin spheres has been truncated ai £ = 6 and the 
parameter Rmt ■ Km ax has been set to 7. The PCPA 
calculations have been performed by a conveniently mod- 
ified version of our KKR-CPA code'^^ aU PCPA cal- 
culations are based on the ASA approximation for the 
site potentials, use several thousands k-points in the full 
Brillouin Zone and 31 energies over a complex integra- 
tion contour. For both LAPW and PCPA calculations, 
the core electrons treatment is fully relativistic while a 
non-relativistic approximation is used for valence states. 
Finally, we present CEF calculationa^i^ with the charge 
multipolar expansion truncated at £ = 0. The concentra- 
tion dependent parameters required by the CEF have 
been obtained from the linear regressions of the 'qV 
data generated from supercells with random occupancies 
and the required mean atomic concentrations and are re- 
ported in Tables |T] and [Tll Depending on which was the 
source of the parameters, the CEF calculations shall be 
referred to as CEF-PCPA or CEF-LAPW. Using the for- 
malism of the previous Section, for both PCPA and CEF 
calculations we set Wi = 1 for all lattice sites, Aqq ~ 1 
and Xim = for £ > 0. 



A. qV laws 

In Sect. Ill we have presented the CEF functional as an 
approximation for the GCPA functional and have shown 
how this is equivalent to assume the linearity of the 'qV 
laws. In Figs. (5) and (6), we plot the 'qV curves from 
our PCPA calculations for the binary bcc Cuo.50Zno.50 
and the quaternary Alo.25Cuo.25Nio.25Zno.25 fee random 
alloys. It is surprising to observe how much accurately 
the PCPA data can be fitted by straight lines. The cor- 
relations coefficients obtained from the linear regression 
of the same data differ from unit by about 10"^. Simi- 
lar very high correlations are always obtained from the 
analysis of PCPA 'qV data, as it is evident by looking 
at Table HI As it is shown in Table HIl also LAPW data 
present high correlations, although the corresponding lin- 
ear fits are not perfect and their correlations deviate from 
unit by 10"^ or 10"'^. This notwithstanding, as argued 
in Sect. Ill C, we believe that the linearity of the 'qV 
relationships within the PCPA is just an approximation. 



TABLE I: CEF parameters obtained by the linear regres- 
sion of the 'qV data from PCPA calculations for random 
CucZni_c alloys in bcc or fee lattices. Ccu and Czn are de- 
fined as the difference between 1 and the correlations obtained 
from the regressions for the Cu and Zn site charges. All the 
quantities are expressed in atomic units. CEF calculations 
using the eoeffieients presented in this Table are referred to 
as CEF-PCPA. 
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TABLE 11: CEF parameters obtained by the linear regres- 
sion of the 'qV data from LAPW calculations for random 
CucZni_c alloys in bee or fee lattices. Ccu and Czn are de- 
fined as the difference between 1 and the correlations obtained 
from the regressions for the Cu and Zn site charges. All the 
quantities are expressed in atomic units. CEF calculations 
using the coefficients presented in this Table are referred to 
as CEF-LAPW. 



e 


acu 


a.zn 


kcu — kzn 


Ccu 


Czn 
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0.445 


4 10"^ 


8 10" 


3 


bcc 
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0.50 
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0.350 


9 10"^ 


3 10" 


3 
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2.317 


0.399 


4 10"^ 
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4 



In order to check out how much accurate it is, we have 
studied one of the most difficult realistic cases, that of a 
high charge transfer ordered alloy, namely the CuZn sys- 
tem. This system has been studied with many different 
theoretical approache o^^'^^i^^'^^i^^'^" . It is also relevant, 
for our present concerns, that the total energy differences 
between fee and bcc geometrical alloy arrangements are 
relatively small. We have executed calculations for all 
the set of 62 bcc and fee based structures reported in 
Refs.lsolandlsil. These structures include several ordered 
crystals for each of the the following Cu atomic concen- 
trations: 0.20, 0.25, 0.33, 0.50, 0.66, 0.75 and 0.80. In 
order to facilitate the comparison, the lattice constants 
have been kept fixed to the values 5.5 and 6.9 a.u., respec- 
tively for bcc and fee based lattices. The results for bcc- 
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and fcc-based alloys are reported in Tables IIIII and IIVI 
respectively. 

The charges from the CEF-PCPA are, in practice, 
identical to those obtained from the PCPA theory for the 
ordered systems. In order to represent the size of these 
tiny differences, we report in Tables [TTTl and ITVl the mean 
square displacement between the two sets of calculated 
charges, < (Ag)^ >. In the worst case, the bcc-based 
DO2 structure, identified in Table [IV] by the number 86, 
we find < (Ag)^ >= 6 lO^'*. Such an excellent agree- 
ment has been obtained for all the set of ordered struc- 
tures considered, in spite of the fact that the CEF input 
has been obtained from random supercells. 

In a previous Letter— we have shown that the CEF 
is able to carefully reproduce the charges from LSMS 
calculations. Moreover, the parameters extracted from 
ordered structure calculations can be used to predict the 
charges for random structures and vice-versa. The qual- 
ity of the CEF predictions was very good either, with 
< (Ag)^ > of the order of 10^^, i.e. about three orders 
of magnitude less than what we have found by the com- 
parison of CEF and PCPA. Since the LSMS calculations 
presented in Ref. 2^ were based on the ASA, we sur- 
mise that the modest lost of accuracy of CEF predictions 
for LSMS wrt. PCPA calculations constitutes a measure 
of the importance of the scattering effects from nearest 
neighbors. These effects, in fact, can be accounted for 
only in a mean field fashion by the PCPA. 

We have also investigated the effects of the spheri- 
cal approximation for the atomic potentials by executing 
full-potential LAPW calculations. In Fig. 8 we plot the 
site charge excesses obtained from LAPW vs. the number 
of unlike nearest neighbors of the same sites, for all the 
structures corresponding to equimolar concentrations. In 
Tables Ull] and [IV] we report the results for < (Ag)^ > 
at all the concentrations. As apparent from Fig. 8, the 
trends of qi are not easily accounted for by the nearest 
neighbors environment only^i, expecially for bcc based 
structures. This notwithstanding, CEF-PCPA calcula- 
tions reasonably account for the LAPW charges. As it 
can be seen in the columns marked as (b) of Tables IIIII 
and IIVI < (Ag)^ > is usually of the order of 10~^, some- 
times less, and about lO""^ in the worst case. In order to 
understand how much these results can be affected by the 
PCPA input coefficients, we have repeated CEF calcula- 
tions by fitting the coefficients from LAPW 'qV data for 
the random alloy configurations corresponding to the rel- 
evant stoichiometrics and reported in Tables IIIII and IIVI 
As shown in the columns (c) of Tables IIIII and IIVI this 
reduces < (Ag)^ > of about one order of magnitude. 

Interestingly, the present CEF-LAPW calculations 
confirm the observations about the transferability of CEF 
parameters in Ref. 2&, where the CEF charges have been 
compared vs. LSMS results. As a typical example, let 
us consider the results for c = 0.50 reported in Table Hill 
The < Aq^ > obtained are alway small: 4 10~^ in the 
worst case and 4 10~^ in the best, corresponding re- 
spectively to structures 71 and 74, while an intermediate 
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FIG. 8: Charge excesses, qi, vs. the number of unlike nearest 
neighbors of the corresponding site, n""' , from LAPW calcu- 
lations for many, bcc and fee based, ordered and disordered, 
configurations, of Cuq. 50 Zno. 50 alloys. Circles and triangles 
represent charges on Cu and Zn sites, respectively. Left frame: 
bcc-based alloys; right frame: fcc-based alloys. The straight 
lines in each panel indicate the best fits obtained by the model 
of Magri et al.— . 
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FIG. 9: Charge excesses on Cu sites qi for the same 
Cuo.5oZno.5o alloys as in Fig. 8. Circles, triangles and squares 
represent the calculated values by CEF-PCPA, CEF-LAPW 
and by the model of Ref. |4l|. In abscissa the charge excesses 
obtained by LAPW calculations are reported. Left frame: 
bcc-based alloys; right frame: fcc-based alloys. In oder to im- 
prove readability we plot also the straight lines q — q^^^^ . 
The deviations from these lines measure the accuracy of the 
various calculations. 



value, 3 10~^ is found for the structure R16, from which 
the CEF-LAPW coefficients have been obtained. The 
same holds for all the concentrations, both for bcc and 
fee structures. A look to the columns (c) in the Tables Hill 
and IIV) in fact, shows that while excellent results have 
been obtained for all the supercell, the random structure 
from which the CEF coefficients have been extracted is 
not necessarily the best performing. 

The above arguments about the charges should not 
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TARLE III: Charge excesses and total energies per atom for bcc-based CucZni-c alloys. The first two columns on the left give, for each system, 
the mean Cu atomic concentration (c) and, when available, the supercell identifier in the database of Refnii (conf). "R" followed by a number, 
e.g., R16, stands for a quasirandom supercell containing the corresponding number of atoms not included in the database. In the third column 
n^jil indicates the mean number of unlike nearest neighbors of Zn sites. The columns from 4th to 7th report the MSD of the charge excesses, 
< (Ag)^ >, obtained by the comparison of different theories: PCPA vs. CEF-PCPA (a), LAPW vs. CEF-PCPA (b), LAPW vs. CEF-LAPW 
(c), LAPW vs. the model of Ref. (d). Columns from 8th to 11th: total energies per atom from PCPA, CEF-PCPA, LAPW and CEF-LAPW 
calculations. The energy zero is given, for each concentration, by the PCPA prediction for the ground state. 
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lead to the conclusion that, for this purpose, any fit is 
comparable with any other. This is clearly shown in Fig. 
9, where the performances of PCPA, CEF-PCPA, CEF- 
LAPW and by the model of Magri et al4i are compared 
for the equiatomic concentration alloys. In the same Fig- 
ure, the distances from the diagonal lines measures the 
differences between the LAPW charges and those by var- 
ious approximations, for all the Cu sites in the supercells 
with c = 0.5. It is there evident that the results by CEF- 
LAPW, marked by open triangles, are much better than 
those by other approximations. 



B. Total energies 

In Tables IIIII and IIVI we compare the total energies 
obtained for CuZn alloys by CEF, PCPA and LAPW 
calculations. We have used the same extended set of bcc 
and fee based structures listed in Ref. [Sll. Since the CEF 
energies contain a, concentration dependent constant, we 



report the quantity AE, defined as the energy difference 
between the structure at hand and the structure that, at 
the same concentration, has the lowest energy, according 
with PCPA calculations. The same AE is plotted in Figs. 
10 and 11 for the Cu concentrations c — 0.25, 0.50, an 
0.75, for which the database of Ref. Isil contains a number 
of structure sufficient to individuate trends. 

Our first observation is that the total energies obtained 
by PCPA and CEF-PCPA calculations perfectly overlap 
on the scale of Figs. 10 and 11, where they are repre- 
sented as filled triangles and open squares, respectively. 
As reported in Tables [1111 and [TVl in fact, the values ob- 
tained by the two methods are different by a few fiRy 
per atom, that is comparable with the accuracy of the 
calculations. Thus, PCPA and CEF-PCPA give indintin- 
guishable results both for the charges (as discussed in the 
previous subsection) and the total energies. Therefore, it 
is compelling to conclude that the CEF theory is a nu- 
merically excellent and powerful tool to reproduce with 
much less efforts GCPA electronic structure calculations. 
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TARLE IV: Charge excesses and total energies per atom for fee-based CucZni-c alloys. The first two eolumns on the left give, for eaeh system, 
the mean Cu atomie eoneentration (c) and, when available, the supereell identifier in the database of Refn-i (conf). "R" followed by a number, 
e.g., R16, stands for a quasirandom supereell eontaining the corresponding number of atoms not included in the database. In the third column 
riunl indicates the mean number of unlike nearest neighbors of Zn sites. The columns from 4th to 7th report the MSD of the charge excesses, 
< (Ag)^ >, obtained by the comparison of different theories: PCPA vs. CEF-PCPA (a), LAPW vs. CEF-PCPA (b), LAPW vs. CEF-LAPW 
(c), LAPW vs. the model of Ref. (d). Columns from 8th to 11th: total energies per atom from PCPA, CEF-PCPA, LAPW and CEF-LAPW 
calculations. The energy zero is given, for each concentration, by the PCPA prediction for the ground state. 
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2.230 


1.986 




16 


10.0 


3 10-'" 


2 10-" 


1 10-" 




0.604 


0.605 


1.054 


0.535 


0.75 


19 


10.0 


5 10-'" 


1 10-* 


5 10-" 




0.729 


0.731 


1.083 


0.709 




22 


12.0 


3 10"" 


2 10-* 


5 10-" 




0.000 


0.000 


0.000 


0.000 




25 


12.0 


2 10"" 


2 10"* 


5 10"" 




0.073 


0.073 


-0.193 


0.057 




27 


6.0 


1 10"** 


5 10"* 


1 10"" 




1.944 


1.949 


3.654 


2.799 




R16 


3.0 


2 10-" 


3 10"* 


2 10-" 




1.088 


1.091 


1.790 


1.409 




30 


10.0 


1 10-" 


6 10"" 






0.544 


0.550 


-4.023 




0.80 


36 


12.0 


2 10-** 


7 10"" 






0.000 


0.000 


0.000 





Moreover, since CEF-PCPA calculations use as an in- 
put the 'qV data obtained from random supercells, the 
perfect agreement obtained for the properties of so many 
different ordered structures that have not been used to fit 
the CEF coefficients has only one possible explanation. 
Accordingly with the discussion in Sects. Ill B and III C, 
both the following conditions must be fulfilled, (i) The 
coherent scattering-path matrix r" of the random alloy 
configuration used as an input must be representative of 
the whole set of ordered structures considered; (ii) the 
linearity of the 'qV laws is almost perfectly observed in 
all the range of values that the charge excesses and the 
Madelung potentials take for the structures considered. 
In Sect. Ill we have offered several arguments support- 
ing the validity of both points above, but we have not 
been able to provide an analytical demonstration. We 
think that the numerical evidence found is very strong 
and compelling. 

Accordingly with the discussion in Sect. Ill, the suc- 
cess of the CEF theory in reproducing the charges or, 
equivalently, the Madelung potentials guarantees the re- 
producibility of any ground state property within the 



GCPA theory through Eq. (|¥T|) . Hence, even spectral 
properties as the DOS or the Bloch Spectral functions, 
although buried, are contained in the CEF functional 
that, if the input parameters are extracted from a GCPA 
theory, inherits all the good and the bad things of same 
GCPA theory 

The comparison with LAPW calculations is more dif- 
ficult, for two different reasons. In first place, these cal- 
culations do not assume mean boundary conditions for 
the wave-functions and use a procedure equivalent to the 
full calculation of the r matrix. In second place, within 
LAPW calculations, the charge multipole summation is 
truncated at some high £ value. With these clarifications, 
the agreement between LAPW and PCPA or CEF-PCPA 
calculations (there is no reason for discussing the last two 
models separately) is quite good. As a general rule, the 
two set of calculations find the same ground states at the 
concentrations considered. In the few exceptions (that 
correspond to the negative figures in the LAPW columns 
of Tables IIIII and IIV|) the disagreement can be explained 
by the fact that the structures indicated as the ground 
state by the two theories are almost degenerate in energy. 
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FIG. 10: Total energy differences with respect to the PCPA 
predicted ground state, AE, for bee based CucZni_c alloys. 
The labels indicated in abscissa identify the various config- 
urations in the database of Ref. [3ll . R stands for structures 
with randomly generated chemical occupations containing 16 
atoms with mean Cu contents c=0.25, 0.5 and 0.75. Open tri- 
angles, open circles, open squares and filled triangles, indicate 
LAPW, CEF-LAPW, CEF-PCPA and PCPA calculations, re- 
spectively. Lines are a guide for the eye. 
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FIG. 11: Total energy differences with respect to the PCPA 
predicted ground state, AE, for fee based CucZni_c alloys. 
The labels indicated in abscissa identify the various config- 
urations in the database of Ref. [sil . R stands for structures 
with randomly generated chemical occupations containing 16 
atoms with mean Cu contents c=0.25, 0.5 and 0.75. Open tri- 
angles, open circles, open squares and filled triangles, indicate 
LAPW, CEF-LAPW, CEF-PCPA and PCPA calculations, re- 
spectively. Lines are a guide for the eye. 



Also the general trends for the total energies are well re- 
produced, as it is visible in Figs. 10 and 11, though the 
PCPA generally underestimate the energy differences. 

Fitting the CEF parameters from the LAPW 'qV laws 
generally improves the agreement. At variance of what 
found for the charges, however, the improvement is quite 
modest. 



In summary: the CEF appears able to perfectly repro- 
duce GCPA calculations for both ordered and disordered 
metallic systems. The reasons why the agreement is so 
excellent are not yet completely understood. Although 
CEF and GCPA theories are both coarse grained versions 
of the DFT, opposite to what numerical results suggest, 
they are not the same theory. In fact, as discussed in 
Sect. Ill, in order to be coincident to the CEF, GCPA 
theories should (i) exactly observe linear 'qV laws and 
(ii) lead to coherent scattering matrices independent on 
the configuration in a fixed concentration ensemble. For 
metallic alloys, these conditions appear plausible and the 
numerical evidence strongly support the view that both 
are nearly satisfied. However we must highlight that the 
condition (i) is not verified for pathologically high values 
of the Madelung field. The comparison vs. LAPW calcu- 
lations suggest that both coarse grained theories, GCPA 
and CEF, are able to reproduce semiquantitatively the 
total energies of the alloy configurations considered. In 
particular, the results by the coarse grained theories are 
strongly correlated with those by LAPW. This fact is 
better elucidated by Figs. 10 and 11, where configura- 
tions belonging to the same fixed concentration ensemble 
are ordered in such a way to have increasing PCPA to- 
tal energies. If the same ordering was not observed by 
some other method for some configuration, this would 
show up as a local minimum in the corresponding curve. 
The most visible of such events, occurs in Fig. 11 for 
c = 0.75, where the curve corresponding to LAPW calcu- 
lations presents a very weak local minimum at the config- 
uration 25. The examination of Figs. 10 and 11 suggests 
that the coarse grained theories are able to give quali- 
tatively correct predictions about ordering for the alloys 
considered, while the fact that they generally underes- 
timate the corresponding energies could imply incorrect 
estimates of the corresponding transition temperatures. 



V. CONCLUSIONS 

We wish to conclude this paper with a summary and 
a few comments. 

We have introduced the class of the GCPA theories, 
that are characterized by (i) a specific ansatz for the ki- 
netic part of the density functional, which is common to 
all CPA-based theories, and (ii) an external model that 
determines the way in which the atomic effective poten- 
tials should be reconstructed and the statistical weights 
to be assigned each. The GCPA class of approximations 
includes most existing CPA-based density functional the- 
ories, to mention a few: the CPA prototype, i.e. the 
single site CPA^i^, the Screened Impurity Model CPA 
(SIM-CPA)ii>i^, the Polymorphous CPA (PCPA)*!, the 
CPA including Local Fields (CPA-hLF)^, the Non Lo- 
cal CPA (NL-CPA)S2. The ansatz (i) consists in apply- 
ing averaged boundary conditions at the surfaces of each 
scattering volume and naturally leads to algorithms re- 
quiring a number of operation that scales as N. As it 
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is discussed by Abrikosov and Johansson^Q, CPA-based 
approximations allow for a careful picture of the spectral 
properties of metallic alloys. The so much criticized re- 
sults of the SS-CPA about the total alloy energies can 
be healed by external models that consider the charge 
distribution in the system. We have shown how this can 
be done systematically by writing the relevant energetic 
contributions as a series involving the charge multipole 
moments in each scattering volume. The truncation er- 
rors of the same series are probably already quite small 
when only the first term is included, as in the case of 
spherical approximations. 

We have derived an expression of the GCPA density 
functional that, together with the above multipole sums, 
includes local 'atomic' terms, completely determined by 
the atomic number of the ion in the volume, and by the 
geometry of the same volume. The local term at the i-th 
site is coupled to the others only through the coherent 
scattering matrix jf and the Madelung potential at the 
same site. Although this kind of coupling, that we have 
called marginal coupling, is not necessarily weak, never- 
theless, it is analytically tractable and it is the source 
of the 0{N) scaling in GCPA theories. We have demon- 
strated that in a GCPA theory all ground state properties 
within a specific sample are functions of the appropri- 
ate coupling Madelung potential only, or, equivalently, 
of the charge multipole moments at each lattice site. To 
put it into other words: we have demonstrated that the 
GCPA approximations realize a coarse graining of the 
Hohenberg-Kohn density functional, since only a part of 
the information conveyed by the electronic density field, 
namely the charge multipole moments, is actually enter- 
ing in the GCPA approximate functional. Moreover we 
have suggested that the explicit form of the GCPA func- 
tional dependence on the multipole moments can be ob- 
tained in a fixed concentration ensemble by the numerical 
integration of the 'qV relationships for a random alloy 
configuration belonging to the same ensemble. The above 
procedure does not rely on the linearity of the 'qV laws. 

We have re-derived the CEF^'' as a sensible approxi- 
mation of the GCPA theories, with which it would co- 
incide provided the 'qV were exactly linear as claimed 
by many groups. The present derivation allows for the 
inclusion of higher order multipole moments. A very re- 
markable feature of the CEF theory is that it shares the 



same structure of the MST. In fact, the minimization of 
the CEF requires the solution of a set of Euler-Lagrange 
equations that has the same structure of the Korringa- 
Kohn-Rostoker (KKR) matrix at zero energy and wave- 
vector^°'^^. More specifically, as it can be seen by com- 
paring Eqs. (HI) and (|52p. the site-diagonal response func- 
tions, fti^ii', and the Madelung coefficients, My^^i', in 
the CEF theory, play the role of the site diagonal scat- 
tering matrices and the KKR structure constants in the 
MST theory. The correspondence is not only formal, 
since the ai^LU are single-site quantities in the same sense 
of the SS scattering matrices^*^ and, in plain analogy with 
them, are related with the SS response to the appropriate 
perturbing field^^. 

In the present paper we have provided several formal 
arguments and strong numerical evidences that CEF and 
GCPA theories lead to very similar results, the discrep- 
ancies being of the order of the numerical errors. We 
have also shown that CEF and GCPA theories are able to 
reproduce the charges and the total energies for many or- 
dered alloy configurations. In our view the coarse grained 
theories, GCPA and CEF, constitute a valuable alterna- 
tive to full DFT calculations. 

Although the CPA theory was proposed many year ago 
with the purpose of dealing with substitutionally disor- 
dered alloys, we think that we have shown that today 
GCPA theories are able to deal with ordered intermetal- 
lic compounds too. Therefore, the fact that CPA-based 
theories are able to cope with sophisticated model of dis- 
order is not an original sin but, rather, an added value. 

The computational performances of the CEF have been 
discussed in more details elsewhere^^. Here we like to 
mention that the possibility of evaluating total energies 
for thousand atoms in a few seconds CPU time could 
constitute a substantial enlargement of the domain of 
the applications of the DFT. 
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Matrices in the angular momentum space are denoted by a 
single underline, double underline is used for matrices both 
in the angular momentum and in the site spaces. Angular 
momentum components are denoted by capital letters, L = 
{£,m), L' = {£',m'), . . . , site components by small Latin 
letters, i,j, 

In the present paper the convention is used that for any 
function of the position fi{r) stands for /(r — R,) with 
r e Vi. 

In our notation {Vi'^^'^} stands for the set of all the 



Madelung coefficients of the i-th site, {V/j^f", V/; 



and {1/^^^^} for the set of all the Madelung coefficients at 
all crystal sites, i.e. {F*^-*^} = {V^^'^jU {Fif -^^l U • • • . 
The notations {qi} and {q} have a similar meaning. 



